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Abstract 

Lord  Rayleigh  derived  the  degree  of  polarization  for  singular  molecular  scattering  phenomena  in  the 
atmosphere  as  a  function  of  the  relative  angular  orientations  of  the  source,  scattering  center  and 
observer.  This  paper  extends  the  Rayleigh  model  by  giving  the  azimuth  angle  for  the  linearly 
polarized  scattered  irradiation  detected  by  the  observer.  Together  these  two  parameters  give  a 
complete  Stokes  vector  characterization  of  Rayleigh  scattering  in  the  earth-horizon  coordinate 
system.  An  experimental  validation  confirms  the  analytical  results.  The  validation  methodology 
uses  a  polarization  axis  finder  attached  to  a  COTS  (commercial-off-the-shelf)  digital  camera.  Each 
image  supplies  enough  information  to  calculate  the  skylight  polarization  azimuth  angle  in  real  time. 
The  equation  for  the  skylight  polarization  azimuth  angle  is  used  to  create  a  pictorial  representation  of 
the  polarized  sky  by  using  colors  to  represent  azimuth  angles.  A  color  scheme  is  used  to  display  the 
degree  of  polarization  and  the  polarization  azimuth  angle  of  a  vehicle  at  several  different  times  in 
daylight.  An  analysis  of  the  reflected  polarized  light  from  the  vehicle  is  included. 
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cos^  = 


sin  or  cos  9-  sin  6  cos  a  cos  (<p  -A) 


Equation  (4)  does  not  include  interactions  between  the  incident  solar  radiation  and  large  atmospheric 
particles  (cloud  droplets,  ice  crystals,  aerosols,  etc.)  and  effects  due  to  multiple  scattering  and 
molecular  anisotropy. 

Validation  Methods 

Figure  4  is  part  of  an  extensive  experimental  validation  study  confirming  the  accuracy  of  Eq.  4.  The 
experimental  apparatus  consists  of  a  polarization  axis  finder7  attached  to  a  digital  camera.  The 
polarization  axis  finder  is  a  linear  polarization  element  with  its  transmission  axes  tangent  to 
concentric  circles  radiating  from  its  center,  as  shown  in  Fig.  5.  It  attenuates  the  radiation  into  two 
wedge-shaped  regions  which  are  aligned  parallel  to  the  polarity  of  the  incident  linearly  polarized 
light.  Figure  4  shows  the  superposition  of  a  circular  protractor  with  its  center  coincident  with  the 
center  of  the  wedge  image.  A  line  through  the  protractor  corresponds  to  the  ip-value  obtained  from 
Eq.  (4)  for  the  appropriate  time  of  day.  This  line  should  pass  through  the  center  of  the  wedge-pattern 
formed  by  the  polarization  axis  finder  confirming  the  theoretical  predictions  of  our  model  relating  to 
Eq.  4.  Several  quantitative  methods  can  be  used  to  validate  Eq.  (4).  One  method  completely 
eliminates  the  need  for  using  the  Stokes  parameters.  Using  Matlab8,  the  x-coordinate  of  the 
minimum  value  along  each  row  (y-coordinate)  of  pixels  can  be  calculated  for  the  wedge  region.  The 
slope  of  the  best-fit  linear  equation  for  these  minimum  values  corresponds  to  the  polarization 
azimuth  angle.  A  more  exact  method  for  determining  the  accuracy  of  Eq.  (4)  uses  the  Stokes 
parameters  calculated  from  the  axis  finder  data. 
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Calculating  the  Stokes  Parameters  from  Axis  Finder  Data. 

The  four  Stokes  parameters,  S0,  Si,  S2  and  S3  for  completely  polarized  light  propagating  along  the 
+z-axis  are9 

So  —  Eox2  +  Eoy2  Si  =  Eox2  —  Eoy2  S2  =  2  Eox  Eoy  cos  8  S3  =  2  Eox  Eoy  sin  8  (5) 

where  Eox  and  Eoy  are  the  instantaneous  amplitudes  of  the  two  orthogonal  components  of  the  electric 
field  vector  and  8  is  their  phase  difference.  The  polarization  azimuth  angle  (ip)  is  determined  from10 

(6) 

It  will  be  shown  shortly  that  three  of  the  Stokes  parameters  (S0,  Si,  S2)  are  required  to  determine  ip, 
even  if  the  light  is  elliptically  polarized. 

If  light  is  transmitted  through  a  linear  polarizer  with  its  transmission  axis  oriented  at  an  angle  0  with 
respect  to  the  x-axis,  the  intensity  of  the  transmitted  light  can  be  expressed  as 

1(0)  =  Vi  (So  +  cos  20  Si  +  sin  20  S2)  (7) 

1(0)  denotes  an  intensity  measurement  corresponding  to  a  particular  value  of  0.  Substituting  the 
angles  0°,  45°  and  90°  for  0  into  Eq.  (7)  confers  that 

So  =  1(0) +  1(90)  S,  =  1(0)  - 1(90)  S2  =  2  1(45)  -  S0  (8) 
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Figure  6  shows  how  Eq.  (8)  can  be  used  in  conjunction  with  a  polarization  axis  finder  to  determine 
the  Stokes  parameters  So,  St  and  S2.  The  transmission  axes  of  the  axis  finder  are  circular  and 
tangents  to  them  give  the  orientation  of  an  infinitesimal  linear  polarizer.  All  points  along  the  line 
between  A  and  A*  in  Fig.  6  correspond  to  linear  polarizers  with  0  =  0°.  All  points  along  the  line 
between  B  and  B*  and  C  and  C*  correspond  to  linear  polarizers  with  0  =  45°  and  0  =  90° 
respectfully.  Figure  7  shows  the  result  of  this  methodology.  The  triangular  markers  on  this  graph 
correspond  to  the  average  skylight  polarization  azimuth  angle  obtained  from  the  axis  finder,  using 
the  Stokes  parameters.  They  are  all  nearly  coincident  with  circular  markers  which  were  obtained 
from  Eq.  (4).  Similar  results  from  a  multitude  of  different  viewing  directions  result  in  the 
conclusion  that  Eq.  (4)  is  an  excellent  tool  in  predicting  skylight  polarization  azimuth  angles. 

The  circular  field  of  view  of  the  axis  finder  was  measured  to  be  10.2°  and  the  corresponding  circular 
image  had  a  radius  of  480  pixels.  However,  a  cropped  circle  of  radius  140  pixels  was  used  to 
calculate  the  skylight  polarization  azimuth  angle.  This  reduction  in  the  original  field  of  view  to  3° 
produces  a  very  narrow  range  of  skylight  vg-values  and  improves  the  accuracy  of  the  validation 
methodology. 

Visualizing  the  Skylight  Azimuth  Angle 

In  previous  work1  *'12,  we  describe  how  the  Stokes  parameters  can  easily  be  encoded  in  a  daylight 
scene.  We  assign  RGB  pixel  values  to  the  normalized  values  of  Si,  S2  and  S3  at  each  pixel  site  in  the 
scene  as  follows: 
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R=int[127.5(l-S0]  G  =  int[  127.5  (1  -  S2)]  B  =  int[127.5  (1  -  S3>] 


(9) 


The  pseudo-color  scheme  described  in  that  work  closely  relates  to  the  Poincare  representation  of 
polarized  light.  The  equator  of  the  Poincare  sphere  plays  a  very  special  role  in  our  pseudo-coloring 
scheme.  It  corresponds  to  linearly  polarized  light  (S3  =  0)  and  is  also  used  to  encode  the  polarization 
azimuth  and  ellipticity  angles  in  a  daylight  scene,  as  given  in  Eq.  (10).  Substituting  x  for  vy  in  Eq. 
(10)  produces  a  color-mapping  scheme  for  the  x-images  as  given  in  Eq.  (1 1). 

R  =  int[  127.5  (1  -  cos  2vp)]  G  =  int[l  27.5  (1  -sin  2vp)]  B  =  127  (10) 

R  =  int[  127.5  (1  -  cos  2X)]  G  =  int[127.5  (1  -  sin  2X)]  B-127  (11) 

The  degree  of  polarization,  P,  varies  between  0  and  1 .  The  simplest  method  of  encoding  this 
parameter  in  a  daylight  scene  is  to  use  the  equation 

[R,G,B]  =  255  P  (12) 

A  pictorial  representation  of  a  semi-hemispherical  region  of  the  celestial  sphere,  as  shown  in  Figures 
8d,  8e,  8f  and  8j,  8k  and  81,  was  created  in  an  Excel  spreadsheet.  Matlab  scripts  that  incorporate 
Eqs.  (2,  3  and  12)  for  the  P-parameter  and  Eqs.  (4  and  10)  for  the  ^/-parameter  were  written  to 
colorize  the  areas  between  the  grid  lines.  The  pseudo-colored  images  relate  to  clear  skies  on  April 
26,  2000  at  north  latitude  42.5  degrees  and  west  longitude  83.0  degrees.  The  skylight  polarization 
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parameters  are  encoded  into  the  drawings  as  they  would  appear  looking  toward  the  earth,  as  from  a 
satellite  above  the  earth. 

The  Azimuth  Angle  in  a  Daylight  Scene 

The  vehicle  images  in  Fig.  8  were  acquired  using  an  Epson  850Z  digital  camera.  In  previous  work12 
we  described  a  methodology  to  obtain  the  four  Stokes  parameters  for  each  scene  pixel  using  a  COTS 
digital  camera.  In  this  paper  we  present  some  applications  of  these  tools  and  techniques  to  the 
analysis  of  the  polarization  azimuth  angle  (vg)  for  daylight  illuminated  scenes.  Since  ip  depends  on 
the  first  three  Stokes  parameters,  the  analysis  is  independent  of  S3  and  consequently  wavelength 
independent.  Many  factors  are  important  to  characterize  polarization  phenomena  for  the  imagery  in 
Fig.  8  including  the  ambient  light  sources,  material  properties  of  the  scene  objects,  and  relative 
orientation  of  the  camera,  the  object  and  the  light  source. 

Daylight  scenes  contain  unpolarized,  broad  band  sunlight,  partially  polarized  skylight,  and  indirect 
diffuse  reflection  from  the  ground.  The  polarization  characteristics  of  these  three  sources  are  quite 
different.  Direct  solar  illumination  is  nearly  collimated  and  originates  from  a  localized  source. 
Skylight  illumination  is  a  distributed  source  where  maximum  P  originates  from  a  great  circle  90° 
from  the  sun.  The  ground  plane  reflection  is  very  unpolarized  and  in  nearly  all  cases  decreases  the 
average  P  for  the  illumination  incident  on  the  vehicle  surface. 

The  optical  axis  of  the  camera  used  to  obtain  Fig.  8  points  toward  the  north,  and  it  is  horizontal  to 
the  ground  plane.  The  normal  to  the  vehicle  vertical  panels  facing  the  camera  points  south.  The 
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vehicle  surface  is  diffuse  and  reflects  incident  light  into  a  very  large  solid  angle  toward  the  camera. 
The  specular  reflection  from  sunlight  off  these  panels  is  primarily  directed  into  the  ground  plane. 
Since  the  sky  is  a  distributed  source,  light  originating  from  the  southern  horizon  can  enter  the  camera 
from  specular  reflection  off  of  the  vertical  panels.  All  other  reflected  illumination  entering  the 
camera  from  the  vehicle  must  be  due  to  diffuse  reflections. 

Equation  13  predicts13,14  that  the  angle  of  incidence  for  sunlight  incident  on  the  vehicle  is  given  by 
the  following  expression: 

cos  i  =  sin  z  sin  A  sin  s  sin  y  +  sin  z  cos  A  sin  s  cos  y  +  cos  z  cos  s  (13) 

where  i  =  angle  of  incidence,  z  =  zenith  angle  of  sun,  A  =  azimuth  angle  of  sun,  s  =  slope  of  plane 
surface  and  y  =  azimuth  angle  of  plane  surface.  The  azimuth  of  the  plane  is  measured  westward 
from  south  and  is  zero  when  the  normal  to  the  plane  points  toward  the  south.  Equation  13  is  derived 
from  the  vector  dot  product  of  a  unit  vector  R  towards  the  sun  and  a  unit  vector  n  normal  to  a  plane 
surface.  Analysis  of  Fig.  9  shows  that 

R  =  sin  z  sin  A  i  +  sin  z  cos  A  j  +  cos  z  k  (14) 

n  =  sinssinyi  +  sinscosyj  +  cossk  (15) 

At  07:00  the  sun  is  near  the  eastern  horizon  with  a  =  3.8  and  A  =  254.9  degrees.  Thus  skylight  is 
highly  polarized  near  the  meridian,  opposite  the  side  of  the  vehicle  facing  the  camera,  which  is 
shown  in  Fig.  8d.  Also,  the  skylight  vy-value  along  the  horizon  is  approximately  90  degrees  (see  Fig. 
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2a).  Since  the  sun  is  located  in  the  northeastern  quadrant  of  the  sky,  no  direct  solar  irradiation  is 
incident  on  the  panels  facing  the  camera  since  the  vehicle  is  positioned  parallel  to  an  east-west 
direction. 

Figures  8a  and  8g  show  that  paint  color  affects  the  polarization  state  of  the  reflected  light.  A  good 
example  of  this  phenomenon  is  the  door  panel,  which  has  two  distinct  paint  colors.  The  reflected 
light  from  those  regions  have  different  values  for  P  and  \\i.  The  ig  angles  are  approximately  60  and 
90  degrees  for  the  two  regions. 

At  19:00  the  sun  is  near  the  western  horizon  in  the  northwestern  quadrant  of  the  sky  with  a  =  14.8 
and  A  =  95.4  degrees.  Similar  to  the  previous  case,  no  direct  solar  radiation  is  incident  on  the 
southern  panels.  Since  the  relative  positions  of  the  sun,  vehicle  and  camera  are  nearly  symmetrical 
with  respect  to  the  N-S  direction  at  19:00  compared  to  07:00,  the  P  and  images  at  these  two  times 
are  very  similar.  Figure  8  also  shows  that  the  azimuth  angle  of  the  reflected  light  from  the  vehicle 
panels  is  approximately  equal  to  the  skylight  azimuth  angles  near  the  meridian  and  the  horizon,  i.e. 
in  the  direction  of  the  digital  camera. 

Figure  10  shows  the  laboratory  experiment  used  to  obtain  the  Stokes  parameters  from  vertical 
panels.  These  studies  indicate  that  diffuse  reflections  originating  from  incident  polarized  light 
produce  azimuth  angles  that  vary  for  different  pigments.  One  of  the  pigments  showed  that  diffuse 
reflections  have  reflection  azimuth  angles  that  are  the  mirror  images  of  the  azimuth  angles  of  the 
incident  light.  This  phenomenon  was  observed  for  the  following  conditions:  (1)  viewing  along  the 
normal  to  the  vertical  panel,  (2)  polarized  incident  light,  and  (3)  all  incident  and  azimuth  angle 
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configurations  for  the  linearly  polarized  light  source.  This  behavior  is  clearly  seen  in  Fig.  1 1  and  the 
daylight  scenes  in  Fig.  8.  This  data  also  shows  that  the  P-value  (Fig.  12)  is  very  dependent  upon  the 
azimuth  angle  of  the  incident  light  for  angles  of  incidence  near  90  degrees. 

An  interesting  relationship  exists  between  the  polarization  azimuth  angles  produced  by  sunlight 
reflected  off  of  vertical  vehicle  panels,  and  the  azimuth  angles  of  the  incident  skylight  oriented  along 
the  meridian  and  on  the  horizon.  These  angles  are  equal  as  shown  in  the  next  paragraph. 


If  r  is  a  unit  vector  1 80°  with  respect  to  an  incident  sun  ray,  then 

r  =  sin  z  sin  A  i  +  sin  z  cos  A  j  +  cos  z  k  (16) 

Let  n  be  a  unit  vector  parallel  to  the  surface  normal  of  a  vehicle  panel.  Then,  n  =  j  and  the  vector 
cross  product  produces 

R  =  r  X  n  =  -  sin  a  i  +  cos  a  sin  A  k  (17) 

The  geometry  in  Fig.  13  shows  that 


cos  (3  = 


sin  a 

-v/l- cos2  a  cos2  A 


(18) 


The  polarization  azimuth  angle  is  obtained  from  \\j  -  180  -  p. 


Using  the  12:00  sun  position  in  Eq.  (18)  produces  a  reflected  solar  ip-value  of  155.7  degrees. 
Inserting  tp  =  0  and  A  =  0  into  Eq.  (4)  yields 


cos  i//  = 


sin  a 

Vl  -  cos2  a  cos2  A 


(19) 
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Equation  (19)  gives  the  skylight  9-value,  as  seen  from  the  earth,  of  a  point  on  the  meridian  (A  =  0) 
and  on  the  horizon  (9  =  0).  Skylight  traveling  from  this  point  toward  a  vertical  panel  on  the  vehicle 
propagates  along  the  normal  to  the  panel.  It  reflects  off  the  panel  as  a  mirror  image  of  the  incident 
skylight.  The  camera  records  a  9-value  that  is  the  supplement  of  the  value  obtained  from  Eq.  (19), 
namely  155.7  degrees  at  12:00.  According  to  the  color  key  for  9-values  in  Fig.  8,  this  9-value 
corresponds  to  a  medium  green  color.  This  is  the  color  seen  in  Fig.  8h  for  vertical  panels  with  a 
surface  normal  toward  the  south.  The  Stokes  parameters  for  reflected  light  from  the  vehicle  panels 
(averaged  over  10  X  10  pixel  areas)  also  gives  9-values  near  150  degrees. 

Temporal  registration  errors,  originating  from  the  30  second  data  acquisition  period  of  the  Epson 
850Z,  are  not  significant  in  obtaining  valid  9  values.  The  change  in  luminance  due  to  solar  radiation 
and  skylight  during  the  30  second  data  acquisition  period  also  has  little  effect  on  the  Stokes 
parameters  for  these  targets. 

The  Degree  of  Polarization  in  a  Daylight  Scene 

Figures  8a  and  8c  at  07:00  and  19:00  show  much  less  contrast  and  higher  P-values  than  the  imagery 
at  12:00.  Figures  8d  and  8f  at  07:00  and  19:00  show  that  highly  polarized  skylight  exists  at  all 
altitudes  in  the  southern  sky,  i.e.,  along  the  meridian  and  directly  opposite  the  vehicle.  Figure  14 
shows  that  both  edges  and  vertical  panels  can  reflect  highly  polarized  light  into  the  camera.  Since  P 
is  large  for  both  of  these  regions,  there  is  little  contrast  between  them  in  the  imagery.  Since  the  sun 
is  behind  the  vertical  plane  at  07:00  and  19:00,  direct  sunlight  can  not  contribute  to  producing 
polarized  reflected  light  from  the  vehicle.  At  12:00,  the  altitude  and  azimuth  of  the  sun  are  55°  and 
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319°  respectfully.  Highly  polarized  skylight,  visible  to  the  panels,  is  near  the  eastern  and  western 
portions  of  the  sky;  hence,  its  contribution  to  the  degree  of  polarization  of  reflected  light  from 
vertical  panels  along  the  optical  axis  of  the  camera  is  much  less  than  at  07:00  and  19:00. 

As  Fig.  15  illustrates,  only  solar  specular  reflections  from  rounded  edges  can  enter  the  camera  at 
12:00  while  specular  reflections  from  vertical  panels  are  directed  into  the  ground  plane  and  not 
visible  to  the  camera.  However,  sunlight  that  reflects  diffusely  from  both  edges  and  vertical  panels 
can  be  seen  by  the  camera.  Laboratory  studies'5  show  that  the  highest  degree  of  polarization  (0.55) 
occurs  for  obtuse  angles  between  the  incident  unpolarized  light  and  the  light  reflected  from  the 
surface.  Since  obtuse  diffuse  reflection  angles  are  directed  into  the  ground  plane,  edges  will  show  a 
larger  degree  of  polarization  than  vertical  panels.  This  phenomena  explains  why  there  is  more 
contrast  between  the  vertical  panels  and  edges  of  the  vehicle  in  the  12:00  P- image. 

Conclusions 

This  paper  gives  a  complete  description  of  linearly  polarized  light  originating  from  single  molecular 
scattering  mechanisms.  Since  there  is  no  elliptical  polarization5  from  primary  Rayleigh  scattering, 
the  Stokes  vector  for  this  case  contains  two  independent  components:  the  degree  of  polarization  first 
derived  by  Lord  Rayleigh,  and  the  azimuth  angle  described  in  this  paper.  Both  of  these  parameters 
are  written  with  respect  to  an  earth-horizon  coordinate  system.  The  latter  is  convenient  for 
predicting  polarization  effects  due  to  solar  illumination  in  planetary  atmospheres.  Both  parameters 
are  dependent  upon  the  relative  orientation  of  the  illumination  source,  scattering  centers  and 
observer.  The  degree  of  polarization  is  important  for  determining  the  presence  of  linearly  polarized 
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light  in  daylight  scenes.  The  azimuth  angle  is  particularly  important  for  determining  the  orientation 
and  curvature  of  secondary  scattering  surfaces  that  are  illuminated  by  direct  sunlight  and  indirect 
polarized  skylight. 

Temporal  registration  errors  caused  by  30  second  time  delays  in  the  image  acquisition  process  do  not 
appear  to  be  important  for  obtaining  accurate  vy- values  in  solar  illuminated  imagery.  The  small 
variations  in  solar  luminance  during  the  data  acquisition  period  also  have  little  effect  on  the  accuracy 
of  Stokes  parameter  measurements  for  these  scenes. 
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Figure  Captions 


Figure  1.  The  observation  plane  PAO  containing  an  incident  solar  ray  PA,  a  scattering  center  P,  the 
observer  O  and  the  scattering  angle  Q. 

Figure  2.  Scattered  waves  that  produce  linearly  polarized  skylight  with  polar  axes  arranged  in 
concentric  circles  about  the  direction  of  propagation. 

Figure  3.  The  geometry  for  finding  the  polarization  azimuth  angle  for  primary  Rayleigh  scattering. 
Figure  4.  Polarization  of  skylight  at  07:31  on  Sept.  7,  2001  at  Payson,  AZ  (long.  =  111 .35°  W,  lat. 
34.23°  N);  a  =  17.3°,  A  =  274.8°,  0  =  20°,  cp  =  0°.  Equation  4  predicts  that  \\i  =  75.2°. 

Figure  5.  A  polarization  axis  finder  is  a  linear  polarization  element  with  its  transmission  axes 
tangent  to  concentric  circles  radiating  from  its  center. 

Figure  6.  The  geometry  of  a  polarization  axis  finder  used  to  obtain  the  Stokes  parameters. 

Figure  7.  Comparison  of  skylight  polarization  azimuth  angles.  Values  marked  with  circles  are 
obtained  from  Eq.  (4).  Experimental  values  that  are  obtained  from  the  axis  finder  are  marked  with 
triangles. 

Figure  8.  The  images  above  relate  to  clear  skies  on  April  26,  2000  at  north  latitude  42.5  degrees  and 
west  longitude  83.0  degrees.  A  normal  to  a  vertical  vehicle  panel,  which  faces  the  digital  camera, 
points  south;  the  camera  view  is  toward  the  north.  The  skylight  drawings  are  pictorial 
representations  of  the  polarized  sky  in  the  southern  semi-hemispherical  region  of  the  celestial  sphere. 
The  skylight  polarization  parameters  are  encoded  into  the  drawings  as  they  would  appear  looking 
toward  the  earth,  as  from  a  satellite  above  the  earth.  The  degree  of  polarization  (P)  is  encoded  into 
the  scenes  using  [pixel  value]  =  255  P.  The  polarization  azimuth  angle  (ip)  is  encoded  into  the 
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scenes  using  RGB  pixel  values  corresponding  to  R  =  int[127.5  (1-cos  2y),  G  =  int[127.5  (1-sin  2\j/) 
and  B  =  127.  Skylight  v(/-values  were  obtained  from  Eq.  (4). 

Figure  9.  Geometry  for  deriving  the  angle  of  solar  incidence.  R  is  a  unit  vector  towards  the  sun  and 
n  is  a  unit  vector  normal  to  a  surface. 

Figure  10.  In-house  experimental  setup  to  obtain  the  Stokes  parameter  for  reflected  light  from  a 
vertical  panel. 

Figure  11.  Reflected  ip- values  obtained  from  the  in-house  experiment  shown  in  Figure  10.  p  is  the 
azimuth  angle  of  the  incident  light. 

Figure  12.  Reflected  P-values  obtained  from  the  in-house  experiment  shown  in  Figure  10.  p  is  the 
azimuth  angle  of  the  incident  light. 

Figure  13.  Geometry  for  finding  \[i.  vp  =  1 80  -  p. 

Figure  14.  At  07:00  and  19:00,  highly  polarized  skylight  along  the  meridian  is  incident  on  the 
target. 

Figure  15.  At  12:00,  only  specular  reflections  of  sunlight  from  edges  can  enter  the  camera;  specular 
reflections  from  vertical  panels  are  toward  the  ground. 
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Figure  1 
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Figure  5 
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Polarization  Azimuth  Angle  of  Skylight 
Viewing  South  at  Altitude  =  20° 
Latitude  =  34.2°  N,  Longitude  =  111.4°  W 


Figure  7 
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(a)  P  of  Vehicle  at  07:00 


(d)  P  of  Skylight  at  07:00 


(g) y  of  Vehicle  at  07:00 


(j)  v|/  of  Skylight  at  07:00 


Color  key  for  v|/-values: 


(b)  P  of  Vehicle  at  12:00 


(c)  P  of  Vehicle  at  19:00 


(e)  P  of  Skylight  at  12:00 


(h)  y  of  Vehicle  at  12:00 


(k)y  of  Skylight  at  12:00 


(f)P  of  Skylight  at  19:00 


(i)  y  of  Vehicle  at  19:00 


(1)  y  of  Skylight  at  19:00 
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Figure  11 
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Degree  of  Polarization  from  a 
Vehicle  Panel  in  a  Vertical  Plane 
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Figure  12 


30 


Figure  13 
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